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We study ground state and finite temperature properties of disordered heavy fermion metals by 
using a generalization of dynamical mean field theory which incorporates Anderson localization 
effects. The emergence of a non-Fermi liquid metallic behavior even at moderate disorder is shown 
to be a universal phenomenon resulting from local density of states fiuctuations. This behavior is 
found to have a character of an electronic Griffiths phase, and can be thought of as a precursor of 
Anderson localization in a strongly correlated host. The temperature dependence of the conducting 
properties of the system reveal a non-trivial interplay between disorder and inelastic processes, which 
are reminiscent of the Mooij correlations observed in many disordered metals. 
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I. INTRODUCTION 

The interplay of disorder and strong correlations re- 
mains one of the least understood topics of contempo- 
rary condensed matter physics. These effects are believed 
to bear relevance to many problems that have attracted 
recent attention, such as the metal-insulator transition 
(MIT) in two-dimensional electron systems (i Disorder ef- 
fects are also likely to be important for the understand- 
ing of the puzzling non-Fermi liquid (NFL) behavior of 
several heavy fermion compounds!^ In some of these sys- 
tems, impurities seem to play only a subsidiary role: the 
explanation for the anomalous behavior is more likely 
to be found in the physics of quantum criticalityj^i^i^iSii 
even though a complete description is still lacking 
However, in other heavy fermion systems, disorder seems 
to play a more essential role and seems to be at the origin 
of the NFL behavior,iSiiiiiii& 

Several attempts have been made to address theoret- 
ically the role of disorder in heavy fermion compounds 
(see an overview below). Many experimental results can 
be described within the so-called Kondo disorder model 
(KDM)^^ or, equivalently, the dynamical mean field the- 
ory (DMFT) of disordered Kondo/ Anderson lattices r^tii 
at least above the lowest temperatures. Essential to 
this description is the consideration of the full distri- 
bution of local Kondo temperatures Tk- For sufficient 
disorder, it has a large weight as Tk — > describing 
the presence of dilute 1ow-Tr- spins that dominate the 
thermodynamic and transport properties. However, the 
KDM/DMFT predictions relied on a fine tuning of the 
bare disorder that suggested that a more accurate mi- 
croscopic foundation was necessary. Conspicuously miss- 
ing in this scheme were fiuctuations in the conduction 
electron local density of states (DOS), a quantity that 
is crucial for the determination of Tk- This was reme- 
died by two of us through a generalization of the DMFT 
that incorporates such Anderson localization effects while 
keeping its local treatment of correlationsjiS^ As a re- 
sult, the question of the extreme sensitivity to the bare 
disorder was solved. One important result of this study 



is the emergence of a power-law distribution of Kondo 
temperatures P {Tk) (x T^~^, where a depends continu- 
ously on the strength of disorder W and decreases as W 
is increased. As the distribution becomes more singular, 
several thermodynamic quantities become divergent, in a 
manner characteristic of Griffiths phasesiSS The system 
eventually localizes at a critical disorder strength Wmit- 
Since there is no magnetic phase transition this has been 
dubbed an electronic Griffiths phase. 

This initial wor k^^i^^ employed the slave boson large- 
N theorjiSiiS^ to solve the auxiliary single-impurity prob- 
lems posed by the method. While versatile, powerful and 
yet computationally cheap, this approach presents some 
disadvantages, the main one being the difficulty of work- 
ing at finite temperatures. It should be reminded that, 
since we deal with wide distributions of Tk, we need to 
be able to describe well the full crossover from T <C Tk to 
T ^ Tk, which is not possible with the slave boson large- 
N treatment. In particular, conspicuously missing are 
inelastic scattering processes. Besides, though giving a 
good description of the low-energy Fermi- liquid regime of 
the single-impurity problem, this treatment does not in- 
corporate high- and intermediate-energy incoherent pro- 
cesses. Another impurity solver is therefore needed to 
assess the importance of these intrinsically finite-T and 
finite-energy features. A particularly useful method, able 
to fill this gap at a reasonable computational cost, is sec- 
ond order perturbation theory in U. We have used this 
method to both complement and cross-check the slave 
boson results. 

Besides the results of the initial workfiSiiS which have 
been confirmed by both methods, some of our main 
conclusions are: (i) there is a subtle interplay between 
conduction and f-electron site disorder that leads to a 
surprising non-monotonic dependence of the conducting 
properties on disorder, a feature that is likely unique to 
Kondo/ Anderson as opposed to Hubbard models; (ii) lo- 
calization effects are essential for the determination of the 
distribution of Kondo temperatures and a KDM/DMFT 
description is clearly insufficient, especially if one starts 
from an experimentally measured discrete distribution; 
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(iii) the interplay between disorder and inelastic pro- 
cesses can lead to a temperature dependence of the con- 
ducting properties that is reminiscent of the ones found 
by Mooij and others in several strongly correlated disor- 
dered metals 

This paper is organized as follows. We review the 
disorder-based mechanisms of non- Fermi liquid behavior 
in the next subsection. Section ^describes the model of 
disordered Anderson lattices we studied and the meth- 
ods we employed to solve it. Section ITTTI focuses on the 
detailed results obtained within the slave boson large-N 
method. This expands considerably on the previously 
published resultsii^iiS In Section ITVl we show the results 
obtained with perturbation theory. Finally, we wrap up 
with a general discussion of the strengths and limitations 
of this study and point out possible future directions in 
SectionEl Some details of the computational procedures 
are given in an Appendix. 



A. Brief overview of disorder-based mechanisms of 
non-Fermi liquid behavior 

1. Kondo-disorder models and the electronic Griffiths phase 

The KDM was proposed early on to account for the 
temperature dependence of the Cu nuclear magnetic res- 
onance (NMR) line widths in VCn^-xP^x {x = Q.^-l)P 
It assumed that disorder in a heavy fermion system gen- 
erates random spatial fluctuations of the exchange cou- 
pling constant J between local moments and conduc- 
tion electrons (the Kondo coupling). Each local moment 
was assumed to undergo the Kondo effect in a manner 
that is completely uncorrelated with the others and each 
with a characteristic energy scale, its Kondo tempera- 
ture Txj- Even narrow Kondo coupling distributions 
lead to a wide distribution of Kondo temperatures due 
to the latter's exponential dependence on the former. As 
a result, at low temperatures, many spins are quenched 
while a few percent remain unquenched and dominate, 
giving rise to singular, NFL thermodynamic properties 
(specific heat and magnetic susceptibility). The NMR 
results in VCvl^-xPAx {x = 0.5 — 1) are well described 
within this picture if the distribution function P (Tk) is 
such that P {Tk 0) -> constii^ The KDM gained a 
natural theoretical setting within the DMFT— of a dis- 
ordered Anderson/Kondo latticeiiS*ii In this approach, 
each conduction electron site exchanges single particle 
excitations with an average "cavity" bath, which is in 
turn self-consistently determined. This treatment be- 
comes exact in the limit of infinite dimensionality and 
is the natural generalization of the Curie- Weiss mean 
field theory of magnets to a fermionic system. Its treat- 
ment of disorder is equivalent to the well-known coher- 
ent potential approximation {CPA)^ In an Anderson 
lattice description, the localized f-electron is hybridized 
with its adjacent conduction electron orbital and spa- 
tial fiuctuations are preserved through the random dis- 



tribution of hybridization strengths. The local moments 
are no longer independent since their distribution self- 
consistently determines the cavity bath. Besides show- 
ing that the KDM corresponds to a rigorous limit of a 
microscopic Hamiltonian, the DMFT enabled the cal- 
culation of other properties such as the resistivit y^^i^T ^ 
dynamic magnetic susceptibilityfiii optical conductivitj*2£ 
and magneto-resistance(S£ with good agreement with ex- 
periments. The non- Fermi liquid behavior of these quan- 
tities hinged on the condition that P [T^ ^ 0) ^ const. 
Subsequent experiments of muon spin rotation^S and 
NMR in high fields^^ showed some inconsistencies with 
the KDM/DMFT, suggesting that inter-site correlations, 
which are absent from that approach, may play a cru- 
cial role at the lowest temperatures. However, annealing 
studies have further emphasized that the consideration 
of disorder effects is indispensable!^ 

Despite its success, the KDM/DMFT description suf- 
fered from a basic deficiency, which can be ascribed to 
its extreme sensitivity to the bare disorder distribution. 
Indeed, the connection between the distribution of J and 
the distribution of Kondo temperatures is too rigid and a 
proper description always relies on fine tuning. In partic- 
ular, discrete distributions of bare parameters can never 
generate a P{Tk) such that P {Tk ^ 0) ^ const.. Like- 
wise, power law distributions of Tk are often necessary 
and it is not clear how they can be obtained within the 
KDM/DMFT. 

More recently, two of us have pointed out that this 
can be solved through the inclusion of localization effects. 
From a basic theoretical point of view, the importance of 
this modification is unquestionable. Disorder scatters the 
conduction electrons giving rise to spatial fiuctuations in 
their wave function amplitude. These Anderson localiza- 
tion precursor effects in turn give rise to fiuctuations of 
the conduction electron local DOS. The Kondo temper- 
atures are exponential functions of the local DOS and 
will show a wide distribution for mild disorder strengths, 
even in the absence of fluctuations in JM^ Besides, direct 
experimental determination of the distribution of J from 
x-ray absorption fine-structure (XAFS) experiments in 
UCu5_j;Pdj; have shown that additional conduction elec- 
tron disorder is necessary for the interpretation of the re- 
sults within the KDM/DMFTi22i2a Finally, the addition 
of localization effects has proved to be just the necessary 
ingredient for the elimination of the extreme sensitivity 
to the bare disorder and for a much more universal de- 
scription. 

The average cavity bath of the DMFT, however, com- 
pletely neglects DOS fiuctuations and a more general 
treatment is necessary. Progress could be made by 
means of the statistical dynamical mean field theory 
(statDMFT), which incorporates the full distribution of 
the conduction electron local DOS, while keeping the 
treatment of local correlations already present in the 
DMFTi2ii5£ The treatment involves solving a fully self- 
consistent loop: the f-electron fiuid gives rise to an effec- 
tive disorder potential for the conduction electrons, while 
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the latter's DOS fluctuations determine the distribution 
of Kondo temperatures. We enumerate the main conclu- 
sions of our analysis tiSii^ 

• Universality: The distributions of several physical 
quantities (Kondo temperatures, local DOS of f- 
and conduction electrons, scattering T-matrices) 
assume a universal log-normal form for weak to 
moderate disorder, irrespective of the form of the 
bare distribution. We have verifled this for Gaus- 
sian, square and discrete bare distributions of con- 
duction electron on-site energies and hybridization 
strengths. This is in contrast to the KDM/DMFT 
results and reflects the mixing of many different 
sites connected by the extended conduction elec- 
tron wave function. 

• Electronic Griffiths phase: Increasing the disorder 
generates wide Tr- distributions, leading to a Grif- 
fiths phase with non-Fermi Hquid behavior. This 
Griffiths phase is not tied to any magnetic phase 
transition but is electronic in origin: it is gener- 
ated by the precursors to the Anderson localization 
transition. 

• Metal-insulator transition: There is an Anderson- 
type localization transition at a critical value of the 
conduction electron diagonal disorder. 

• Non-monotonic conductivity as a function of the 
disorder: The typical conduction electron DOS, 
which vanishes at the localization transition and 
serves as a measure of the conducting properties 
shows a counter-intuitive non-monotonic behavior 
as a function of disorder for a wide range of fillings. 
This surprising feature originates in the interplay 
between the bare and the f-electron disorder po- 
tentials and is tied to the proximity to the Kondo 
insulating (pseudo-)gap. 

2. Magnetic Griffiths phase, scenario 

An alternative theoretical scenario for disorder- 
induced non-Fermi liquid behavior is the magnetic Grif- 
fiths phase In the vicinity of magnetic phase tran- 
sitions disorder fiuctuations induce rare regions with an 
enhanced local critical temperature. These large clusters 
are ordered on the scale of the correlation length and 
act as effective spins. Though rare in occurrence they 
carry a considerable amount of magnetic entropy and the 
overall effect is the appearance of singular, thermody- 
namic responses. This magnetic Griffiths phase picture 
has been advocated as a source of NFL behavior in dis- 
ordered heavy fermion systems 4^ However, very recent 
results seem to point to several difficulties encountered 
when this scenario is applied to experimental systems, as 
follows. 

The entropy problem. The amount of magnetic en- 
tropy observed experimentally in most disordered heavy 



fermion systems seems much too high to be compatible 
with the magnetic Griffiths phase picture. Taking the 
measured specific heat of, say, UCu5_a;Pda;, we estimate 
that about 5% of the sample would have to participate 
in the spin-^ clusters. This implies an average cluster 
separation of 2 — 3 lattice constants, ruHng out cluster 
sizes exceeding this distance. These small clusters sug- 
gest instead that the "unquenched" localized moments of 
the KDM/DMFT or the electronic Griffiths phase offer 
a much more natural explanation. 

Effects of dissipation. Other important limitations of 
the magnetic Griffiths phase scheme have also been em- 
phasized in recent work by Millis, Morr, and Schmalian^S^ 
who have carefully examined the effects of dissipation 
caused by the metallic bath. This work suggests that the 
dissipation caused by itinerant electrons is so pronounced 
that quantum tunneling of even moderately sized mag- 
netic clusters will be suppressed. Although the emer- 
gence of a magnetic Griffiths phase is a well estabHshed 
phenomenon in disordered insulating magnets, this result 
seems to bring into question its relevance to itinerant sys- 
tems. 



3. Spin glass precursors 

Finally, another possibility is to invoke the proximity 
to a spin glass quantum phase transition. Several theo- 
retical schemes predicting non-Fermi liquid behavior in 
the vicinity of a spin-glass quantum critical point have 
been proposedi^Si^i*^ Spin glass phases have been iden- 
tified in the phase diagram of some heavy fermion al- 
loys (UCu5_j;Pdj;, for X > 1.5)^ and structurally disor- 
dered compounds (URh2Ge2)4i More interestingly, ev- 
idence of glassy dynamics in the absence of freezing at 
very low temperatures has been seen in UCu5_a;Pda; 
(x = 1, 1.5)^ and Ce(Ruo.5Rho.5)2Si2r^ with confiict- 
ing results pointing to a very low freezing temperature in 
UCua.sPdi.si^ These experiments seem to suggest that, 
if there is a true spin glass transition, freezing temper- 
atures are strongly suppressed in a wide portion of the 
phase diagram. 



II. THE MODEL AND ITS SOLUTION 

A. The statistical dynamical mean field theory 

A simplified Hamiltonian capable of capturing the es- 
sential physics of disordered metals with locaHzed mo- 
ments is provided by a disordered Anderson lattice 

H^H, + Hf + Hhyb. (1) 
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where 



in imaginary time as 



Hf 

Hhyb 



E + c^E /Ir/^T/li/ji; (3) 

JO- j 

E(^^-4'^^-+H-c.). (4) 



In Eqs. H2I4|I . Cj^ (/jtr) annihilates a conduction (f-) elec- 
tron on site j with spin projection cr, t is the nearest 
neighbor hopping amplitude, ^ is the chemical potential, 
U is the f-site Coulomb repulsion, Ef is the f-energy level, 
and we introduce random conduction electron on-site en- 
ergies (ej) and hybridization matrix elements Vj. These 
are chosen from given distributions Pi (e) and P2 {V) , 
taken to be either square or Gaussian, with width and 
standard deviation W , respectively. We have also stud- 
ied discrete cases of P2 {V). There is a large degree of 
uncertainty as to a realistic model of disorder for heavy 
fermion alloys. A rather thorough study of the local f- 
site environment in the alloys UCuR -.rPdx {x = \ and 
0.5) was carried out in Refs. I32l3.^ through XAFS ex- 
periments. These authors were able to determine the 
amount of Pd/Cu site interchange as well as the U-Cu 
bond length distributions. In order to accommodate both 
types of fluctuations one must, in principle, allow for 
a distribution of both hybridization strengths and on- 
site conduction electron energies. On the other hand, 
when the local moments are randomly replaced by non- 
magnetic elements (the so-called "Kondo holes"), spatial 
fluctuations of Ef should also be includedjii Throughout 
the paper, we use the half bandwidth D as energy unit. 
For the Bethe lattice (to be introduced later), D = 2\f2t 
if the coordination number is 3. This should be con- 
trasted to Ref. [T3 . where W is measured in units of t. 
Note also that the f-level energy Ef is always measured 
relative to the chemical potential. This is to ensure that, 
by making it negative and large enough in absolute value, 
we always work in the local moment/Kondo regime. 



We worked within the framework of the 
statDMFT ]ff"^iffi? This treatment is able to incorpo- 
rate both strong local correlations and Anderson 
localization effects in a fully self-consistent fashion. 
Although the method has been described before in 
the context of the disordered Hubbard modelf2i we 
will briefly review it with the dual goal of setting the 
notation and extending it to the disordered Anderson 
lattice. It starts by focusing on a generic unit cell j of the 
lattice, containing an f-site and its adjoining conduction 
electron Wannier state, and writing its effective action 



Seff (j) = Sc (j) + Sf {j) + Shyb (j) , 

Scij) = E / / dr'cl{r)[S{r^r') 
„ Jo Jo 



(5) 



Sf (j) 



X {dr + ej - m) + Acj (t - r')] c,. (r') , (6) 
''^dr Y.^dr + Ef)fl{T)f,„{T) 



+ Unfj^(T)nfji(T) 



(7) 



Shyb (j) = E / h4 (^) (^) + H- c] , (8) 
^ Jo ^ J 

where nfj^ — f]afj<y In writing Eqs. H5I8|I . a simplifi- 
cation has been made of retaining only quadratic contri- 
butions in fermionic fields after integrating out the other 
sites (besides the instantaneous Hubbard term), much 
like the usual dynamical mean field theoryi^ The bath 
(or "cavity") function Acj (r) in Eq. © is given by 



A,, (r) ^t^Y. ir) , 



(9) 



Lm — l 



where the sum extends over the z nearest neighbors and 



T 



(j) 



is the Green's function for propagation from nearest 
neighbor site I to nearest neighbor site to, calculated with 
the site j removed. Integrating out the remaining con- 
duction electron Cj^, we get the effective action of an 
auxiliary single-impurity Anderson model at each site j 

g. Jo JO 

iSr+Ef)+Af,iT-T')]f,, (r') 

+ / drUuf,^ (r) nf,^ (r) , (10) 
Jo 

where the hybridization function for the f-site is, in Mat- 
subara frequency space. 



Afj {iuj„) 



•^j 



(11) 



The solution of this impurity problem is the major dif- 
ficulty in this treatment. We implemented two different 
methods of solution, which will be expanded upon in the 
next subsections. The aim is to calculate the local f- 
electron Green's function 

G%Hr)^-(T[f,Ar)fl{0)]), (12) 
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under the dynamics dictated by lllOII . It is conveniently 
parameterized by its self-energy 



1 



(13) 



It is also convenient to define a local conduction electron 
Green's function 



(r) 



T 



such that 



cj. (r) c , (0) 



1 



lUJn - ej + H - Acj (lUJn ) - "Pj [li^n ) 



where 



iuJn - Ef - Y.fj (iujn) ' 



(14) 



, (15) 



(16) 



We note that the (iujn) function describes the local 
scattering of the conduction electrons off the f-shell at 
site j, incorporating information about both elastic and 
inelastic processes. 

All the information a generic site j has about the rest 
of the lattice is encoded in the bath function which 
should be viewed as a functional of the conduction elec- 
tron lattice Green's function. Since a fully analytical 
treatment is impossible, we have to solve the equations 
numerically. For this purpose, we formulated the prob- 
lem on a Bethe lattice, where things are considerably 
simplified as explained in Ref. |^ In this case, nearest 
neighbors I and m become disconnected once j is removed 
and only local Green's functions survive 



G[^Li^^n)^6l, 



= Sl,r..G'°'^'^ (ZC.„) . 



Finally, these last objects can be computed from an ac- 
tion at site I in almost all aspects identical to H5I8|I . the 
only difference now being that the bath function sum 
runs over the z — 1 nearest neighbors (here labeled by m) 
only 



= f 



z-1 

E 

m— 1 



(17) 



Note that, on the right-hand side, we do not need to 
specify that site j has been excluded as the removal of 
I completely disconnects sites labeled by m from site j 
(this property is specific to the Bethe lattice). The reap- 
pearance of g1°^(') ( iiUn), whose distribution is identical 
to that of G^;'^^"'^ {iwn) since all sites are equivalent, closes 
the loop and estabHshes a recursive set of stochastic equa- 
tions. When the interaction is turned off {U = 0), this 
treatment reduces to the well-known self-consistent the- 
ory of localizatiouj'*? here generaHzed to a two-band lat- 
tice. When we take the coordination to infinity 2; ^ oo 
keeping t = tj^fz = const., our treatment reduces to the 



DMFT of correlations and disorder In the latter 
case, the disorder treatment is equivalent to the CPAjSi 
which has no localization transition. 

A full solution of Eqs. 1151 17|l for given distributions 
P\ (e) and/or Pi (F) involves solving an ensemble of im- 
purity problems self-consistently. Physically, the conduc- 
tion electrons propagate through a disordered lattice and 
scatter off conduction site potential fiuctuations as well as 
f-site resonances. These resonances, in turn, describe the 
formation of localized moments, whose local Kondo tem- 
peratures fiuctuate as well, refiecting a disordered con- 
duction sea environment. Complete statistical informa- 
tion can in principle be obtained from the distributions of 
the various renormaHzed local quantities. We stress that 
a random distribution of any bare parameter causes all 
renormalized quantities to fiuctuate as a result of the self- 
consistent nature of our treatment. Therefore, even if we 
include only fiuctuations in the conduction sea through 
P\ (e), a distribution of Kondo temperatures ensues^^^ 
This is easily seen from the approximate formula for the 
Kondo temperature in the Kondo Hmit (X^tr "/i^ ~ 1)^ 



Txj ~ D exp 



1 



P3J3J ' 



(18) 



where pj is the local density of states (DOS) seen by the 
f-site 



Pj 



lm[Af,{0-irj)] 



(19) 



and Jj is the Kondo coupling constant, given in the 
Kondo limit by 



1 



1 



\Ef\ \Ef + U\ 



(20) 



Even if Vj is not random, the local DOS is, because of the 
denominator in Eq. lfTT|l . As a result of the strong expo- 
nential dependence in Ijl8ll , even mild localization effects 
can be strongly enhanced and should be seriously consid- 
ered, specially in disordered heavy fermion systems. 



B. The impurity solvers 

An important part of the method we employed is the 
solution of the impurity problems posed by the ensemble 
of effective actions given by Eq. Ijl()|l and its counter- 
part for a site with one nearest neighbor removed. We 
concentrated mostly on two methods of solution, which 
we now briefiy describe: the slave boson, large-N based, 
mean field theory and second order perturbation theory 
in U. In order to unclutter the notation, we drop in 
the next subsections the site index j and the superscript 
loc. Details of the numerical treatment are given in the 
Appendix. 
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1. Slave boson mean field theory 

This method gives a good description of the low tem- 
perature, Fermi Uquid regime of the Anderson impurity 
problem in the limit U ^ oo and is extensively covered 
in the literature iSiiS^ Its main advantage is the ability 
to capture the zero temperature fixed point correctly as 
well as the exponential nature of the low energy scale. 
Its treatment of the self-energy, however, does not incor- 
porate inelastic processes to leading order. Besides, it 
has a spurious phase transition at a finite temperature, 
where in reality there should be only a smooth crossover. 
For these reasons, we confine it to the zero temperature 
limit, where it is a useful guide. As appHed to our prob- 
lem, the method has been described in Appendix D of 
Ref. and we will merely state the results, generalized 
to the Matsubara frequency axis and at T = 0. The local 
f-electron Green's function is given by 



iuj — ef — qAf (iuj) 



(21) 
(22) 



where the last equality defines the local f-electron quasi- 
particle Green's function and the variational parameters 
£/ (renormalized f-energy) and q (quasi-particle residue) 
are determined from the solution of the set of equations 



duj 



Af{iuj)Gf{iuj) = 0, (23) 



/QO 1 
— e'-^Gf (ic^) = 1. 
-OO ^ 



(24) 



Using 



^e-Tm 

TT 



Eqs. II23I24|I simplify to 

duj 



2 / —Re 

/O 



q + 2 — Re 



Af{iLu)GY{iuj) = Ef-ef, (25) 
duj , 



Gf (ILU) 



= 0. 



Eq. H16|l becomes in this approximation 

qV^ 



$ (iuj) 



luj-ef 



(26) 



(27) 



2. Second order perturbation theory 

The perturbative solution of the single-impurity An- 
derson model with particle-hole symmetry was thor- 
oughly analyzed by Yamada and Yosida3^ii^2iS^ The 
series expansion in U for physical quantities such as the 



specific heat and the spin susceptibility converges very 
fast and even second order results can be usefuli^ Ex- 
tension of the perturbative treatment to the case without 
particle-hole symmetry poses considerable difficulties. A 
particularly useful proposal is the use of an interpolative 
self-energy which recovers the atomic {V 0) and high 
frequency limits Further improvements of the method 
were later suggested i^SiS 

The procedure consists in defining an unperturbed f- 
electron Green's function 



1 



- /i — Af{iLu„ 



(28) 



with a new parameter Jl to be determined later, which 
vanishes at particle-hole symmetry. The interacting 
Green's function is given in II13|I . The interpolative self- 



energy igSSiit 



E/ ztj„ = Un + . s — 

1 — BL,y^>(lU}n) 



where 
and 



T^e-"''G/(ic^„), 



(29) 



(30) 



Gf\r) Gy'i-r). (31) 



(0), 



The last equation l(3T|l is the usual second order diagram 
using the unperturbed Green's function lf28ll for the in- 
ternal lines. The parameters A and B are determined by 
imposing the high frequency and atomic limits, respec- 
tively, and are given bji^ 



A 
B 



n{\ — n) 
no(l - no) 
(1 -r^)^7^ 



Ef+fL 



where 



710 = T^e*""''G^ 



(0) 



(32) 
(33) 

(34) 



Different schemes have been proposed in order to fix the 
free parameter At zero temperature, one can 

ensure that the low energy Fermi Hquid behavior is ob- 
tained by imposing the Friedel sum rulei^ This proce- 
dure cannot be easily generalized to finite temperatures, 
however. One option is to fix /I at its zero temperature 
value even at finite temperatures. Alternatively, one can 
require at any temperature^ 



"0, 



(35) 



which makes A = 1. Finally, a third possibility is impos- 
ing Jl = These three alternatives have been rather 
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carefully compared in Ref. at T = and checked 
against exact diagonalization. The first two methods 
were shown to be almost equivalent whereas the third 
one is inferior. Moreover, comparisons at finite temper- 
atures with Quantum Monte Carlo results confirmed the 
adequacy of imposing Eq. Specific appHcations to 

a clean Anderson lattice model further corroborated this 
conclusion Thus, our results were based on imposing 
condition (jSU- It should be remembered, however, that 
the perturbative solution predicts a characteristic energy 
scale that is quantitatively incorrect at large values of 
U, since it is unable to capture the correct exponential 
dependence. Nevertheless, for moderate interactions, it 
still gives reasonable results. Within its Hmitations, this 
perturbative scheme is a relatively fiexible low-cost tool 
to tackle the impurity problem with the great advantage 
of being able to naturally account for inelastic processes. 

We note that a direct comparison between the slave 
boson mean field theory results and second order pertur- 
bation theory is not possible because the former is limited 
to the U ^ oo limit, which is obviously outside the region 
of validity of the latter. The main interest of an analysis 
of both methods, however, resides in the exploration of 
the importance of inelastic processes, which are absent in 
the slave boson mean field treatment. 



III. SLAVE BOSON MEAN FIELD THEORY 
RESULTS 



an order parameter for localization 

A useful measure of this mobility is given by the typ- 
ical value of the local escape rate. This is encoded 
in the imaginary part of the local conduction electron 
Green's function (the local DOS) at zero frequency, pcj = 
ifm [G^°'^ (0 — iS)] . We will, from now on, drop the su- 
perscript denoting the removal of a nearest neighbor so 
as to lighten the notation. As shown originally by An- 
derson, the typical value of the local DOS vanishes when 
the electrons are locaHzed and can be viewed as an order 
parameter for the localization transitioni^ A convenient 
way of accessing the typical value is furnished by the ge- 
ometric average 

p*''^ = exp{h^}, (36) 

where the overbar denotes a disorder average. By con- 
trast, the arithmetic average 

pT^M (37) 

is finite at the transition. A thorough analysis of the 
critical behavior of the local DOS distribution in the non- 
interacting Bethe lattice localization problem was carried 
out in Ref. 

In Fig. n] we show the typical conduction electron DOS 
as a function of disorder for several values of the chemical 
potential. 



We now present the results obtained at T = using the 
slave boson mean field theory as an impurity solver. Most 
of our results were obtained for a uniform distribution of 
on-site conduction electron energies 



In Section rrflCl we also show results for a discrete distri- 
bution of hybridization strengths P2 {V). 



A. Conduction electron typical density of states 

To understand the overall behavior as a function of dis- 
order, it is instructive to consider the transport proper- 
ties of the conduction electrons. Since there are no inter- 
actions among them in our model, their behavior is that 
of a disordered non-interacting electron system. There 
are two sources of disorder, as can be seen in Eq. ifTfill : 
fiuctuations of the local on-site energies Cj and of the f- 
shell resonances described by $j {i(^n)- They are not in- 
dependent, however, since they are inextricably tied by 
self-consistency. Their combined effect acts to decrease 
the conduction electron mobility. 



1. Typical density of states: 



2. Proximity to the Kondo insulator 

We can identify three qualitatively distinct 
behaviorspi^ For /.t < —0.4 and for /i > 0.5, p'^^ is 
a monotonically decreasing function of disorder. For 
/i ~ 0.1, the clean system is a Kondo insulatoi42iSi and 
Pc" = P^c^ = 0- disorder is introduced in the Kondo 
insulator, p'^^ initially increases, reaching a maximum at 
about W « 1.5, after which it decreases monotonically. 
Finally, for -0.3 < ^JL < Q.l and 0.1 < < 0.4 p^yp 
initially decreases, passes through a minimum around 
W ~ 0.5—1, then increases up to a maximum at 
about W ~ 1.5 and eventually becomes monotonically 
decreasing. For all values of /u, the typical DOS vanishes 
at a disorder-induced metal insulator transition (MIT) 
at Wmit ~ 4.5. 

These distinct behaviors can be traced back to how 
close the clean system is to the Kondo insulator poin*£2iSi 
at p ~ 0.1.1- jf TffQ start from the clean insulator, the 
introduction of disorder acts to create states inside the 
gap, thus increasing the DOS at the chemical potential. 
This increase continues until the gap is essentially washed 
out and the system becomes a bad metal. After that, 
localization intervenes and p^^P starts to decrease towards 
the MIT. 

For fillings close to but not at the Kondo insulator 
point, the clean system is a heavy fermion metal at T = 0. 
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Figure 1: Typical conduction electron density of states as a 
function of disorder strength for several values of the chem- 
ical potential /i, using slave boson mean field theory as the 
impurity solver. We used Ef = —1, V = 0.5, except for the 
thick solid line in (a), which is for V = 0. 



The f-resonances ($ (w)) coherently scatter the conduc- 
tion electrons creating a strongly renormalized Fermi liq- 
uid. In the slave boson treatment, <I> (w) diverges at 
Lu = ef, see Eq. lf27jl . corresponding to the limit of uni- 
tary scattering, with a maximally allowed phase shift 
S = 7r/2. We can view its value at the Fermi level 
<i>(0) = ~qV'^/ef (which is real at T = 0) as an effec- 
tive potential coming from the f-electrons. The closer the 
system is to the Kondo insulator, the larger the value of 
4> (0) , the insulator being signaled by the divergence of 
this quantity (or equivalently by e/ — 0). The effect of 
disorder is to immediately start generating spatial fluc- 
tuations of the f-resonances, with different phase shift 
values at the chemical potential. Proximity to the in- 
sulator impHes large, random, almost unitary scattering 
potentials. As a result, metallic coherence is efliciently 
destroyed and the typical conduction electron DOS is 
strongly suppressedpi£*i^ The important role played by 
the unitary scatterers was emphasized in Refs. [13 and 
where the distribution P (1/ [<i> (0)]) was directly com- 
puted and its weight at 1/$ (0) = was shown to corre- 
late with the destruction of coherence. 

There is another equivalent way of understanding these 



Figure 2: Typical conduction electron density of states as a 
function of disorder strength for different values of hybridiza- 
tion V, using slave boson mean field theory as the impurity 
solver {Ef = —1 and /i = 0). 

effects. For small dopings away from the Kondo insula- 
tor, carriers are introduced at the edges of the valence 
or conduction bands deflned by the Kondo insulator gap, 
which have a small DOS (in the Bethe lattice, band edges 
have a square root shape as in three dimensions). As 
has been known for a long time, a region of small DOS 
is particularly sensitive to localization effects introduced 
by disorder. 

As in the previous case, further increase of disorder 
acts to wash out the nearby Kondo pseudo-gap and the 
behavior then becomes very similar to the disordered 
Kondo insulator. We thus have a region with a rather 
non-intuitive increasing p**'^, which can be ascribed to 
the proximity to the Kondo insulator fixed point. The 
behavior at fillings well away from the Kondo insulator 
is much less influenced by the pseudo-gap, see Fig. ^ for 
/i = ±0.5. Although there is a rapid initial decrease of 
p^J^P, followed by a much slower dependence, the typical 
DOS does not exhibit the unconventional increase with 
disorder observed at other flllings. 

3. Role of the hybridization strength 

It is interesting to note that the critical value of dis- 
order for the MIT Wmit depends on the hybridization 
strength. In Fig. [2 we show the disorder dependence 
of the typical conduction electron DOS for different val- 
ues of V. There is hardly any change in Wmit as we 
go from V = to V = 0.5 (cf. also Fig. QJ. However, 
for V = 0.75, the critical disorder strength is clearly en- 
hanced. This figure also illustrates the non-trivial nature 
of the self-consistency. Indeed, the two types of disor- 
der coming from fluctuations in ej and $j are clearly not 
independent, since the addition of f-site disorder as we 
turn on V from to 0.75 acts to increase the mobility for 
W > 2.8. The self-consistently determined solutions of 
the impurity problems effectively help screen the conduc- 
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Figure 3: Distributions of the logarithm of Im [Gc"" (0 — iS)] 
for several values of disorder strength, using slave boson mean 
field theory as the impurity solver: (a) /x = —0.5, (b) = 
and (c) n = 0.1 (V = 0.5, Ef = -1). 



tion electron disorder. Note also how an increased value 
of V pushes the "dip-hump" structure to higher values of 
disorder. Since the Kondo insulating gap increases with 
the hybridization strength, this is consistent with our ex- 
planation for the nature of this non-monotonic behavior. 



Trpc = Im [G'-°'^ (0 — iS)] . We show this by plotting 
P (log Im [Gjf^ (0 — id)] ) for several disorder strengths 
and different chemical potential values in Fig. El (we 
use "log" for the base 10 logarithm). It follows from 
the definition of Trpl^^, Eq. . that it is obtained by 
raising 10 to the power of the average of this distri- 
bution (we use powers of 10 for ease of computation). 
For weak disorder in the metallic cases (Figs. and 
b), P (log Im [G^""^ (0 — iS)] ) is approximately Gaussian, 
even though the bare disorder is uniform, a feature shared 
by several physical quantities. ^^-^^ This is due to the pres- 
ence of correlations between many distant lattice sites 
mediated by the extended conduction electron wave func- 
tion, which introduces a sort of averaging effect. In 
the Kondo insulator case (Fig. Et), however, the distri- 
bution is not Gaussian at weak disorder. Keeping in 
mind that Pc = (log Im [G^"" (0 - i6)] -oo) in the 
clean Kondo insulator, it is clear that the introduction of 
weak disorder has to generate weight at very small irpc- 
Indeed, P (log Im [G^""^ (0 — iS)] ) shows a divergence at 
Im [G^°= (0 - iS)] « 10-3 for W = 0.35. 

For large values of W the distribution becomes ex- 
tremely broad, spanning many orders of magnitude. In 
the case p = —0.5 (Fig. Et), corresponding to a system 
well away from the Kondo insulator filling, the distribu- 
tion broadens and its maximum steadily shifts towards 
lower values as disorder is increased. This is to be ex- 
pected from the monotonic behavior of p*^^. Likewise, at 
p = 0, the non-monotonic behavior of the typical value 
is also clearly reflected in P (log Im [G^°^ (0 — iS)] ) (see 
Fig.Et and compare it to Fig.QJi). 

As we saw, at the Kondo insulating chemical potential 
p = 0.1 and for W = 0.35, the distribution shows a diver- 
gence at Im [G^°^ (0 — iS)] « 10"'^. A similar diverging 
tendency is observed at p = (Fig. Eb) and W = 0.5. 
This is precisely the disorder value where the minimum 
of /O**'^ occurs (cf. Fig. QJi) and which we have been as- 
cribing to the presence of many unitary scatterers due 
to the nearby Kondo insulator. The similarity between 
the two distributions strengthens further our case for the 
importance of the proximity to the Kondo insulator. Ad- 
ditionally and consistent with this, the divergence is to- 
tally absent at p = —0.5, where the role played by the 
Kondo insulator fixed point is much less important. 

It is also interesting to observe in Fig. Efc, how the 
Kondo gap is washed out by disorder: at W = 1.4, where 
pI^p peaks (Fig.^), most of the weight of the distribu- 
tion is already at sizeable values of the DOS and its shape 
is very similar to the metallic cases. 



B. Distribution of the conduction electron local 
density of states 

One of the great advantages of the present approach 
is the possibility of monitoring complete distribution 
functions. Many of the features exhibited in Figs. 
and 12 can be directly read off the distribution of 



C. Distribution of Kondo temperatures 

We now proceed to the analysis of the physical prop- 
erties related to the ensemble of impurity problems. As 
shown beforeiSiiS the distribution of Kondo temperatures 
of the various f-sites is log-normal for weak disorder, but 
broadens and acquires a power-law shape at intermediate 
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Figure 4: (a) Distribution of the logarithm of the Kondo tem- 
perature for a discrete bare distribution of hybridizations, 
taken from Ref. 32. using slave boson mean field theory as 
the impurity solver (here, we use t = 0.5, Ef = —5.54 and 
fj, = —0.2); (b) Comparison between the smooth distribution 
of Kondo temperatures (dot-dashed line) obtained in the stat- 
DMFT and the discrete results of DMFT (vertical solid lines) 
and the Kondo disorder model (vertical dashed lines). 



values ofW~ 0.35 — 0.7. Once this power-law becomes 
singular enough, a Griffiths phase is entered with diverg- 
ing thermodynamic responsesiiS*iS 



1. Universality at weak disorder 

We have noticed that for weak disorder, the shape 
of the distribution of various quantities, including the 
Kondo temperature, is universal, irrespective of the 
shape of the bare distribution of disorder. 

A nice illustration of this effect is given by the case 
where the bare disorder is a discrete distribution. As an 
example, we take the discrete distribution of hybridiza- 
tion strengths, P2 {V), determined in Ref. fsl from XAFS 
measurements in UCus-jjPda;. The resulting distribu- 
tion of logTx is shown in Fig. 0^. It is continuous and 
has a log-normal shape. In Fig. |^ we show the distri- 
bution of Tj( (smooth dot-dashed fine) and compare it 
to the discrete distribution obtained in the DMFT (ver- 



Figure 5: Distribution of the logarithm of the Kondo temper- 
ature for different values of disorder, using slave boson mean 
field theory as the impurity solver {V — 0.5, Ef = —1 and 
fi = 0.2). 

tical solid lines), which corresponds to the limit of infi- 
nite coordination. We also include in the figure the re- 
sults of the Kondo disorder modeli^i^ (KDM) (vertical 
dashed lines), which is very similar to the DMFT. The 
only differences between the KDM and the DMFT are 
that, in the former, no self-consistency is imposed and 
a Kondo instead of an Anderson lattice model is used. 
The difference between the results of the KDM/DMFT 
and the statDMFT is striking. The fiuctuations of the 
conduction electron wave functions incorporated in the 
StatDMFT smooth out the discrete results of the DMFT 
into a universal continuous form. A description of the 
NFL behavior within the KDM/DMFT theory would be 
clearly impossible. This comparison also shows that this 
level of hybridization disorder alone is not able to gener- 
ate non-Fermi liquid behavior even in statDMFT, since 
the distribution of Kondo temperatures goe s to ze ro as 
Tk ^ 0, a point that was stressed in Refs. How- 
ever, if disorder in the conduction electron sites, Pi (e), 
is also included a singular behavior can be obtained (not 
shown) . The inclusion of conduction electron disorder is 
reasonable in VCu^-xPdx, since the Cu-Pd interchange 
affects both Vj and ej. 

2. Emergence of the Electronic Griffiths phase 

In order to identify the emergence of the Griffiths 
phase, we next study the evolution of the distribution 
of Kondo temperatures as the width W of Pi (e) is var- 
ied (with no disorder in V). Typical results are shown 
in Fig. El corresponding to ^ = 0.2. As the disorder 
increases, we find that the overall width, but most signif- 
icantly, the size of the low-Tjf tail rapidly grows. These 
tails assume a power-law form P{Tk) ~ T^~^ , with the 
power a{W) being a monotonically decreasing function 
of disorder. Once again, this behavior cannot be obtained 
in the rigid scheme of the KDM/DMFT without unjus- 
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Figure 6: Exponent of the power-law distribution of Kondo 
temperatures as a function of disorder for different values 
of the chemical potential {V — 0.5, Ef — —1). The hori- 
zontal dashed line indicates the critical value for the emer- 
gence of NFL behavior, where P (Tk) oc const, and x{T) oc 
C(T) /Toe In (To/T). 



tified fine tuning. The thermodynamic response of the 
system assumes a singular, non-Fermi Hquid form as soon 
as a < 1, which happens for sufficiently strong disorder 
W > Wnfl ~ 0.35 — 0.45. Since this behavior does not 
reflect any thermodynamic phase transition, it assumes 
the character of an electronic Griffiths phase. Here, sin- 
gular behavior emerges due to the presence of exponen- 
tially rare events (in our case Kondo spins) vi^hich nev- 
ertheless provide an exponentially large contribution to 
thermodynamic and transport properties and thus dom- 
inate the macroscopic behavior of the system. The W- 
dependence of the exponent a can be easily obtained by 
fitting the tails of these distributions; a representative 
behavior for several values of the chemical potential is 
shown in Fig. El 



3. NFL and the proximity to the Kondo insulator 

In the previous section we have seen how the proximity 
of the Kondo insulator plays a crucial role in determining 
the disorder dependence of the locaHzation and transport 
properties, and leads to the surprising "bad metal" behav- 
ior for a wide parameter range. This was argued to reflect 
the enhanced sensitivity to disorder of those electronic 
states that are very close to the Kondo insulator, and 
which are most easily affected by localization effects. We 
have also established that NFL behavior also emerges as 
a result of disorder- induced density of states fluctuations. 
It is then natural to ask how sensitive this emergence of 
NFL behavior is to the proximity to the Kondo insula- 
tor, which in the clean limit emerges only in a narrow 
parameter range, close (in our case) to ^ 0.1. 

To address this question, we have systematically in- 
vestigated the evolution of P{Tk) as a function of the 
distance to the Kondo insulator, i.e. as a function of the 



Figure 7: Distribution of the logarithm of the Kondo temper- 
ature for different values of the chemical potential and fixed 
disorder W — 0.18, using slave boson mean field theory as the 
impurity solver {V — 0.5, Ef = —1). 



chemical potential ^. The behavior for weak disorder 
(W = 0.18) is shown in Fig.H Despite what one would 
naively expect, these result clearly demonstrate that the 
distributions are the broadest far from the Kondo insula- 
tor. As we can see on this flgure, the distributions narrow 
down as the Kondo insulator is approached from either 
side. As a result, we may expect that the critical disor- 
der strength Wnfl necessary for the emergence of NFL 
behavior should increase closer to the Kondo insulator. 
This surprising result is conflrmed by examining the ^- 
dependence of the exponent a as shown in Fig. El As we 
can see there, for a given W , the exponent a is indeed 
smaller, and Wnfl decreases for larger /i (far from the 
Kondo insulator). 

At flrst sight, these flndings seem in contradiction to 
what one may expect, since we have found that the typ- 
ical density of states decreases close to the Kondo insu- 
lator. Naively, one could then expect the Kondo tem- 
peratures to be depressed as well, leading to broader dis- 
tributions and enhanced NFL behavior. On the other 
hand, we know that the Kondo temperature remains fi- 
nite within the Kondo insulator, despite the fact that 
the density of states at the Fermi energy vanishes there. 
Although surprising at flrst sight, this curious feature 
of Kondo insulators is at present well understood. It 
is called the "strong coupling Kondo effect"iS^iSSi££ It re- 
flects the fact that the Kondo screening is not determined 
only by electronic states precisely at the Fermi energy, 
but also by all the states in an energy interval of order Tk 
around the Fermi energy. Indeed, the average value of the 
Kondo temperature (see Fig. [jjl is the highest precisely 
near the Kondo insulator. In the presence of disorder, 
the value of the Kondo temperature is determined by a 
certain weighted average of the density of states over this 
extended energy interval. When localization is present, 
only the states closest to the Kondo gap band edge will 
be appreciably affected, but since not only those states 
determine Tk, the net effect is washed away. We thus 
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Figure 8: Distribution of the logarithm of the Kondo temper- 
ature for different values of the chemical potential and fixed 
disorder W — 1.4, using slave boson mean field theory as the 
impurity solver iy — 0.5, Ej = — 1). 



Figure 9: Distribution of Kondo temperatures on a log-log 
scale, for different values of the chemical potential and fixed 
disorder W — 1.4, using slave boson mean field theory as the 
impurity solver iy — 0.5, Ej = —1). 



conclude that the proximity to the Kondo insulator, in 
contrast to transport, does not have appreciable effect 
on the emergence of the NFL behavior. Indeed, the crit- 
ical value of disorder Wnfl required for the emergence 
of NFL behavior is found to have a remarkably weak /i 
dependence. 

J^. Universality at strong disorder 

At strong disorder we expect the density of states fluc- 
tuations to completely wash out any trace of the Kondo 
gap, and in addition to broaden the conduction band, 
making it very flat and featureless. As a result, all quan- 
tities are expected to have an extremely weak /i depen- 
dence, leading to a more universal behavior of all quan- 
tities. Such behavior is indeed seen at sufficiently strong 
disorder, where the typical DOS curves (Fig.nj are seen 
to merge around W ~ 1.4. Similar behavior is seen 
in Fig. IHl which shows -P(log Tk) for different /i's at 
W = 1.4. This universal behavior is even more striking 
if this distribution is plotted on a log- log scale (Fig. EJ, 
where an almost perfect power-law tail (a ^ 0.2) is seen 
to span several decades for all the values of /i. This is a re- 
markable example of universality generated by DOS ffuc- 
tuations, completely absent in the KDM/DMFT treat- 
ment. 



D. Distribution of the hybridization function 

A key input to the determination of the Kondo temper- 
ature is the hybridization function A/ (icun) of Eq. ifTTll . 
We show in Fig.llOh the distribution of its imaginary part 
calculated at the chemical potential (zu;„ O — iS). Note 
that, for a featureless bath, it appears in the exponen- 
tial of the Kondo temperature formula, Eq. ifTsll . which 
is thus very sensitive to it. It can be seen in Fig.llOb that 



its distribution is very regular and retains a bell-shaped 
structure for any W. It also inherits the non-monotonic 
behavior observed in the conduction electron local DOS 
(cf. Fig.ElD). 

It is tempting to try to calculate the distribu- 
tion of Kondo temperatures from the distribution of 
Im [A/ (0 — iS)], by naively applying the Kondo temper- 
ature formula, Eq. ltT8|l . This procedure fails, however: 
the actual distribution of Tk has much lower weight at 
small Tk values than is predicted by the Kondo tem- 
perature formula. The explanation for this failure lies 
in the fact, already alluded to before, that Tk is deter- 
mined by a weighted average of Im [Af {uj — iS)] over a 
region around the Fermi level. This can be glimpsed from 
the strong frequency dependence of the typical (geomet- 
ric average) hybridization function close to the chemical 
potential, as shown, for frequencies on the real axis, in 
Fig.llOb. It shows a robust and well-defined pseudo-gap, 
inherited from the nearby Kondo insulator, and a tiny 
narrow peak at the chemical potential. This narrow peak 
is easy to understand: spatial fluctuations due to disor- 
der give rise to narrow peaks within the pseudo-gap, most 
typically at the chemical potential. However, as we have 
remarked before, the Kondo temperature can be flnite 
even if the density of states is zero or almost zero at the 
chemical potentialm^'^^^ In this case, the spectral weight 
right at the chemical potential is unimportant for the de- 
termination of the Kondo temperature. It is dominated 
by a whole range of spectral density away from the Fermi 
level. Since far from the pseudo-gap region the density 
of states is much larger and hence much less affected by 
the spatial fluctuations, the distribution of Kondo tem- 
peratures is narrower than one might guess based on the 
distribution of Im [A/ (0 — iS)] and the Kondo tempera- 
ture formula. 

When considered together, the results of Section IIIII 
show the importance of a self-consistent solution of the 
problem, with a non-trivial interplay between spatial 
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Figure 10: (a) Distribution of the logarithm of Im[A/(0)] 
for different values of disorder and (b) typical value of 
Im [Af {u> — iS)] as a function of frequency a.t W = 3.5, both 
using slave boson mean field theory as the impurity solver. 
The inset in (b) details the behavior close to the chemical 
potential {V = 0.5, Ef = -1 and n = 0.2). 



fluctuations due to localization and strong correlation ef- 
fects. However, an important feature that is missed in 
the slave boson treatment of the impurity problems is the 
presence of inelastic scattering. This will be considered 
in the next section, where we show the results obtained 
with perturbation theory in the interaction. 



IV. PERTURBATION THEORY RESULTS 

We consider now the results obtained at finite T using 
second order perturbation theory as the impurity solver. 



A. Results for a fixed temperature 

Fig. shows the results for the conduction electron 
typical DOS near the Fermi level as a function of the 
disorder parameter W for different values of the inter- 
action energy U and the hybridization V at T = 0.003. 
In order to understand them, we looked at scatter plots 
of realizations of the local effective f-shell potential (real 



Figure 11: Typical density of states near the Fermi surface as 
a function of disorder for different values of the hybridization 
V and the interaction U, using perturbation theory as the 
impurity solver. Other parameters used were T = 0.003, Ef = 
— 1 and 1-1 = 0. 



and imaginary parts of « 0)) and the correspond- 

ing bare potential Cj at the same site, both seen by the 
conduction electrons. These results are presented for two 
different values of the hybridization, V — 0.3 and 0.5, and 
two different values of the disorder parameter, W — 0.7 
and 7, in Figs. El to ^1 The other parameters used in 
the calculation were U = 4, T = 0.003, Ef = —1 and 
/i = 0. The figures present the results obtained in the 
statDMFT calculation (represented by dots) as well as 
the results of DMFT (full line), in which the disorder 
treatment reduces to CPAw^^ It is important to notice 
that we have Im [$(cj « 0)] 7^ (as Im [S/(t^ « 0)] ^ 
at finite T), implying the presence of inelastic scatter- 
ing, a feature absent in the slave boson treatment of 
the last section. Besides, the imaginary part of the self- 
energy gets folded into the real part of $ (w) as well. 
This is a peculiar feature of a two-band model, where 
the effective conduction electron self-energy represented 
by $ (oj) has a real part for which inelastic processes 
contribute. Thus, we should keep in mind that, even 
though in an effective description of conduction electron 
processes Re [$j (w « 0)] and ej are associated with elas- 
tic scattering while Im (w w 0)] is related to inelastic 
processes, both Re[$j(w w 0)] and Im[<I>j((jj w 0)] con- 
tain information on f-electron collisions. 

Let us compare the results for W — 0.7 for both 
V = 0.3 and V = 0.5, which are in Figs. [Hand [H The 
first observation we make is that for V = 0.3 the values of 
Re [$(w ~ 0)] are mainly concentrated around zero, while 
for V = 0.5 they are distributed in a wider range of val- 
ues. On the other hand, the results for Ini[<l)(ci; w 0)] 
show that the inelastic scattering is stronger for V — 0.3 
than for V — 0.5. Concerning the results for V — 0.5, 
the great concentration of sites with large Re [^{uj w 0)] 
explains the great decrease in the typical DOS for low 
disorder shown in Fig.^] These sites act as almost uni- 
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Figure 12: Scatter plot of the effective f-shell potential 
$j(tj ~ 0) and the bare potential ej at the same site (dots) 
for V = 0.3 and W = 0.7, using perturbation theory as the 
impurity solver: (a) real and (b) imaginary parts. The solid 
line is the DMFT result for W = 0.7. Other parameters used 
were U ^ 4, Ef ^ ~1, T ^ 0.003 and ^ 0. 



Figure 13: Scatter plot of the effective f-shell potential 
$j(Lj ~ 0) and the bare potential €j at the same site (dots) 
for V = 0.5 and W = 0.7, using perturbation theory as the 
impurity solver: (a) real and (b) imaginary parts. The solid 
line is the DMFT result for W = 0.7. Other parameters as in 
Fig.im 



tary scatterers, which give rise to a maximally allowed 
scattering phase shift {S = tt/2) for the conduction elec- 
trons, and represent droplets of Kondo insulator within 
the metal, in close parallel to the slave boson results dis- 
cussed of Section Un l ^^i^^ For V = 0.3, although the in- 
elastic scattering is stronger, it cannot compensate for 
the much narrower distribution of Re [$(cj « 0)]. Thus, 
the DOS does not decrease as fast as for V = 0.5. 

As the disorder increases to intermediate values (W ~ 
1.8), the distribution of e becomes larger, causing a 
steady decrease in the typical DOS for V = 0.3. How- 
ever, for V = 0.5, the typical DOS presents a non- 
monotonic behavior similar to the slave boson results for 
—0.3 ^ < 0.1 and 0.1 < ^ < 0.4. As we explained in 
Section Hill this is a result of the fact that, as the dis- 
order increases, the concentration of unitary scatterers 
first increases, then saturates and the bare disorder dom- 
inates over the f- related one ($j {uj))M-^^ Indeed, even 
for V — 0.3 we notice the presence of a slight "shoulder" 
in the typical DOS around W = 1.1, whose origin is the 
same as that of the non-monotonic behavior at ^ = 0.5. 
We call attention to the slave boson results of Fig. |2l 
which similarly show that a decrease of the hybridization 



strength moves the "dip-hump" feature to smaller disor- 
der values. Moreover, as the current calculation was done 
at finite temperature, inelastic scattering also plays some 
role in smoothing out this feature. 

As the disorder continues to increase and the distri- 
bution of e becomes broader, the typical DOS for both 
V — 0.3 and V = 0.5 decreases. These results tend to the 
non-interacting one, as is expected if only the bare disor- 
der plays a role. Indeed, comparing the results for = 7 
(Figs.[niand ll5|l we notice that the realizations for which 
the bare disorder e is large have Re [$(a; « 0)] around 
zero, meaning that the real part of the f-shell disorder is 
not important. Besides, the values of Im [$(w ~ 0)] are 
small for these realizations. In the case of the realizations 
for which e is small. Re « 0)] attains values that are 
larger than the results for W = 0.7 (cf. Figs. HH and [T3|l . 

The above discussion has focused on the conduction 
electron viewpoint. Let us now consider how the pres- 
ence of disorder in e is seen by the f-electrons. For the 
auxiliary one-impurity problem, the important scale is 
the Kondo temperature T^, which measures the coupHng 
between the impurity and the conduction electron bath. 
The presence of disorder in e generates a distribution of 
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Figure 14: Scatter plot of the effective f-shell potential 
$j(aj ~ 0) and the bare potential ej at the same site (dots) 
for V = 0.3 and W — 7, using perturbation theory as the 
impurity solver: (a) real and (b) imaginary parts. The solid 
line is the DMFT result for W = 7. Other parameters as in 
Fig. [El 



Figure 15: Scatter plot of the effective f-shell potential 
^j{uj ~ 0) and the bare potential ej at the same site (dots) 
for V — 0.5 and W — 7, using perturbation theory as the 
impurity solver: (a) real and (b) imaginary parts. The solid 
line is the DMFT result for W = 7. Other parameters as in 
Fig. [El 



Tk, as we have seen before. Thus at finite temperature 
some of the sites have Tk > T, forming a singlet state 
with the bath, while the sites for which Tk < T repre- 
sent an almost free spin, scattering conduction electrons 
incoherently. This incoherence is characterized by a large 
amount of inelastic scattering, signaled by a significant 
imaginary part of the self-energy. Going back to the re- 
sults for W — 0.7, the fact that the inelastic scattering is 
stronger for V = 0.3 than for V = 0.5 refiects the larger 
number of incoherent sites (with Tk < T) in the for- 
mer, since the smaller the hybridization, the smaller the 
Kondo temperature (see Eqs. El and 1201. On the other 
hand, the sites with large Re [$(tj « 0)] for V — 0.5, 
which are responsible for the great decrease in the typi- 
cal DOS, represent sites with Tk > T. 

Fig. El also shows the results for the typical DOS for 
U = 2, and V = 0.5. In this case, the system has particle- 
hole symmetry in the clean limit, presenting a gap in its 
DOS (the Kondo insulator). This explains the fact that 
for small disorder the typical DOS decreases as W de- 
creases. The overall behavior here bears strong similarity 
with the slave boson results of Fig. QJi (circles) and the 
explanation for it has been given in Section Hill 



B. Temperature dependence 

Fig. El presents the results for the typical DOS as a 
function of W for different temperatures. The other pa- 
rameters used were U = A, V = 0.3, Ef = —1 and /i = 0. 
Here is where the interplay between inelastic and elastic 
processes proves to be fairly non-trivial and the use of a 
technique that incorporates both, such as perturbation 
theory, is crucial. First we note that for the lowest dis- 
order value {W — 0.25), the typical DOS decreases with 
increasing temperature. This is made more clear in the 
inset. On the other hand, for W — 1.4, this tendency 
is reversed. Finally, in between these two extremes, the 
temperature dependence can be non-monotonic. A bet- 
ter sense of the overall behavior can be grasped by plot- 
ting the inverse of the typical DOS as a function of T 
for different values of W, as shown in Fig. El It is clear 
that for 0.7 <W< 1.2, the inverse typical DOS shows a 
peak as a function of T, which gradually moves to zero 
temperature as disorder is increased. 

The temperature dependence of the typical DOS can 
be rationalized by looking at the corresponding changes 
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Figure 16: Typical density of states near the Fermi surface 
as a function of disorder for different values of temperature, 
using perturbation theory as the impurity solver. Other pa- 
rameters used were U = 4, V = 0.3, Ef — —1 and /i = 0. 
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Figure 18: Scatter plot of the effective f-shell potential 
$j(Lj ~ 0) and the bare potential €j at the same site (dots) 
for T = 0.0005 and W = 0.7, using perturbation theory as 
the impurity solver: (a) real and (b) imaginary parts. The 
solid line is the DMFT result for W = 0.7. Other parameters 
as in Fig. El 



Figure 17: Inverse typical density of states near the Fermi 
surface as a function of temperature for different values of 
disorder, using perturbation theory as the impurity solver. 
Other parameters as in Fig. Ilfil 



in the distribution of $ (w « 0). For this purpose, we fo- 
cus on the disorder value of = 0.7, for which l/7rp*^P 
has a clear maximum at T « 0.0005. We show in Figs.1181 
and El the scatter plots of $ (w « 0) at T = 0.0001 and 
T = 0.0005, respectively, which should be compared to 
the higher temperature results (T = 0.003) of Fig. El 
From Fig. El we see that the inverse typical DOS in- 
creases as we go from T — 0.0001 to T = 0.0005 and then 
decreases as the temperature is varied up to T = 0.003. 
As can be seen from the figures, this non-monotonic be- 
havior is governed by the effective f-disorder encoded 
in the distribution of Re [$ (w r:! 0)]. Indeed, its vari- 
ance increases from T = 0.0001 to T = 0.0005 but de- 
creases from T = 0.0005 to T = 0.003. Note that the 
imaginary part of <& (a; « 0) always increases with in- 
creasing temperature, reflecting the enhancement of in- 
elastic processes. In terms of T^, this is the same as 
saying that as the temperature increases the number of 



sites with Tk < T, which have a stronger inelastic scat- 
tering, becomes larger. However, in the interval from 
about T — 0.0005 up to T = 0.003, the increase of 
Im [$ {uj ~ 0)] is outweighed by the narrowing of the dis- 
tribution of Re [$ (cj « 0)], which is the dominant contri- 
bution. The above analysis can be similarly extended to 
the other values of disorder shown in Fig. El 

All along we have been using the typical DOS as a mea- 
sure of the conducting properties of the system. This 
is justifled by its interpretation as a escape rate from 
a lattice site and the fact that it vanishes at the lo- 
calization transitiont^i Ideally, one would like to calcu- 
late the conductivity instead. This is a difficult task 
in the present scheme, however, although an approxi- 
mate calculation can be performed, which becomes ac- 
curate close to the localization transitioniSSiSS It requires 
the calculation of the propagator between two different 
sites, which goes beyond our current method, whose fo- 
cus is on local Green's functions only. Even in view of 
all these caveats, however, it is tempting to use the in- 
verse typical DOS as an approximate measure of the re- 
sistivity, specially far from the weakly disordered region 
{W > 0.35). If this is done, then the fanning out of the 
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Figure 19: Scatter plot of the effective f-shell potential 
$j(a; ~ 0) and the bare potential ej at the same site (dots) 
for T = 0.0001 and W = 0.7, using perturbation theory as 
the impurity solver: (a) real and (b) imaginary parts. The 
solid line is the DMFT result for W = 0.7. Other parameters 
as in Fig. 



"resistivity" curves of Fig. El as T — > is reminiscent 
of the Mooij correlations;^ originally observed in disor- 
dered transition metal alloys, but which are also seen 
in heavy fermion alloys, W doped Kondo insulator gT^iT^iTi? 
and even in two-dimensional systems as in the metal- 
oxide-semiconductor field-effect transistors pi Our results 
and the discussion above show that the interplay be- 
tween locaHzation effects and electron-electron interac- 
tions can give rise to the Mooij correlations, without the 
need to invoke other sources, such as electron-phonon 
interactions. ^^•'^4 In addition, we have pointed out how 
the rapid drop in the typical DOS is a consequence of 
the proximity to the Kondo insulator, a region where lo- 
calization effects are particularly large. Taken together, 
these observations point to a close connection between 
Mooij correlations and localization effects in the vicinity 
of a Kondo or a Mott insulator. It would be interesting to 
test these ideas by a direct calculation of the resistivity 
within a scheme capable of incorporating both locaHza- 
tion and strong correlations. 

Finally, regarding the DMFT calculation presented in 
all the figures discussed in this Section, we note that the 
results obtained for Re [$(aj « 0)] fall approximately in 



the region where the concentration of realizations in the 
statDMFT calculation is largest. There is no reason to 
expect DMFT and the statDMFT to give similar results, 
except at weak disorder. Nevertheless, our results show 
that, surprisingly, DMFT can serve as a rough guide for 
the most probable values of Re « 0)]. There is a 
sizeable discrepancy, however, between the DMFT and 
the StatDMFT results for Im [$ (w w 0)] at the highest 
temperature (T = 0.003) and low disorder. The DMFT 
line in this case is an underestimate of the reaHzations 
obtained within the statDMFT. As the temperature is 
lowered a better agreement is obtained. Thus, besides its 
inherent neglect of localization effects, DMFT should be 
used only as a lower bound when gauging the importance 
of inelastic processes in disordered Anderson lattices. 



V. DISCUSSION AND CONCLUSIONS 

We have in this paper extensively characterized the 
physics of disordered Anderson lattices within the stat- 
DMFT scheme, which is able to incorporate both lo- 
caHzation effects and the local correlations coming from 
electron-electron interactions. This was done using both 
large-N methods and perturbation theory for the auxil- 
iary single impurity problems. These are in a sense com- 
plementary approaches. On the one hand, large-N theory 
is ideal for ground state properties, where inelastic effects 
are absent. In particular, it affords a quick and reliable 
way of calculating Kondo temperatures (with the correct 
exponential dependence) and scattering phase shifts at 
T = 0. Without 1/A^ corrections, however, it is unsuit- 
able for a finite temperature calculation. On the other 
hand, second order perturbation theory has the advan- 
tage of being equally flexible with the added bonus of 
incorporating inelastic processes and the temperature de- 
pendence of the scattering phase shifts. Nevertheless, it 
fails to capture the exponential nature of the low tem- 
perature scale of the single impurity problem. Taken to- 
gether, the two methods have enabled us to put on flrm 
grounds the conclusions laid out in previous work )^^i^i? 
namely, the emergence of an electronic Griffiths phase 
in Anderson lattices governed by the proximity to the 
disorder-induced localization transition. In particular, 
several inadequacies of the early Kondo Disorder Model 
(and its formulation as a Dynamical Mean Field The- 
ory) have been given a better theoretical basis once lo- 
calization effects were included. The self-averaging effect 
introduced by the spatial fluctuations of the conduction 
electron wave functions induce a much higher degree of 
universality than is possible in the rigid KDM descrip- 
tion. Furthermore, the perturbation theory treatment 
has also suggested a mechanism behind the ubiquitous 
observation of Mooij correlations in the resistivity of dis- 
ordered materials. 

Even within the conflnes of the approximations of the 
StatDMFT scheme, there are still outstanding issues that 
we would like to resolve. Even though a fully analytical 
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treatment is probably impossible, it might be feasible 
to devise an approximate parameterization of the stat- 
DMFT on the Bethe lattice, specially at weak disorder. 
In this respect, the universality of the distributions of 
"dressed" quantities, such as the various local Green's 
functions and the Kondo temperatures, which are ei- 
ther Gaussian or log-normal is a useful guide. A "toy" 
model that assumes a simple form for the distribution of 
Acj {itOn) can be written down, which recovers the Grif- 
fiths singularities obtained in the numerical treatment 
This "toy" model may prove useful for a calculation of 
the resistivity and for generalizations of the statDMFT 
treatment. 

On the other hand, a specially important effect ne- 
glected in the statDMFT is two-particle inter-site corre- 
lations, particularly in the spin channel. Indeed, the pro- 
liferation of poorly quenched Iow-Tr- spins in our treat- 
ment generates a large amount of entropy that must be 
relieved at low temperatures through inter-site correla- 
tions. Indeed, experimental evidence in favor of spin- 
glass dynamics at low temperatures in UCu5_a;Pda;!^ and 
URh2Ge2^ makes the inclusion of inter-site correlations 
more pressing. A promising avenue of attack would be 
to remain true to the spirit of the DMFT and use its 
extended version. Several treatments along these lines 
have been attempted4i*^Si2& The challenge in our case 
is to incorporate both the dynamical inter-site correla- 
tions of the latter treatments and the spatial fluctuations 
of Kondo temperatures of the electronic Griffiths phase 
we find in our approach. We defer the discussion of this 
problem to a future publication. 
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Appendix A: BRIEF DESCRIPTION OF THE 
NUMERICAL METHOD 

The set of stochastic equations defined in H5I17|I must 
be solved in two steps. First, action (I5I8II is solved with 
the bath function with z — 1 nearest neighbors defined 
in itTTjl . This determines self-consistently the distribu- 
tion of local conduction electron Green's functions with 
one nearest neighbor removed G^"'^^"''' (ia;„). Next, the 
same action (I5I8II is solved, this time with the bath func- 
tion in constructed from the previously determined 
^c°'^'"'' (i^n)- This step involves no self-consistency and 
yields the distribution of G^°'' (»w„) as output. Since the 
latter bath function is a sum over z nearest neighbors its 



statistical fiuctuations are reduced compared to the for- 
mer one. Thus, G^°'^ (ia;„) is more narrowly distributed 

than G^°i^^^'^ («t^n)- Yet, we expect the qualitative behav- 
ior to be the same. We have therefore focused on the first 
step of the procedure only. 

For a given impurity solver, this disordered Bethe lat- 
tice problem was solved for z = 3, by sampling the dis- 
tribution of G^°'^'''' (iw„) from an ensemble of N sites, as 
proposed originally in Ref. We have generally used 

= 70 — 100 since we have checked that the results do 
not change by taking TV = 200. We have thus determined 
the distribution of various local properties. 

The equations were solved on a discrete mesh along 
the Matsubara axis. The mesh is set by the Matsubara 
frequencies in the perturbative treatment at finite tem- 
peratures and by an arbitrary finite discrete mesh in the 
slave boson mean field theory at T = (up to 32000 
points). The choice of the imaginary frequency axis is 
due to a greater numerical stability. When the disorder 
is strong, the various Green's functions show large fiuc- 
tuations. However, these are much more pronounced on 
the real frequency axis, where they give rise to several 
peaks and gaps. 

The slave boson treatment consists in finding the two 
mean field parameters q and £/ by solving the set of two 
non-linear Eqs. H25I26II . For that, we used the Powell 
hybrid method. The integrals were calculated with stan- 
dard adaptive quadrature routines. Since we used a fi- 
nite frequency mesh, it was important to extrapolate the 
value of GY (*t^) with the asymptotic form 1/iuj for val- 
ues of w greater than the largest mesh value. In this fash- 
ion, a wide range of Kondo temperatures is covered, going 
down to almost machine precision in the clean metallic 
case. When disorder is present a given impurity prob- 
lem may not have a solution even at T = 0, because the 
strong spatial fiuctuations may cause the local DOS to 
vanish at the Fermi level. In this case, there are two pos- 
sible regimes for the impurity, which have been carefully 
analyzed iS^iSSiSi The analysis shows that there is a criti- 
cal coupling constant Vc such that the ground state is a 
singlet for V > Vc (the so-called "strong coupling Kondo 
effect"), whereas the local moment remains unquenched 
\{ V < Vc- When a solution could not be found, this 
corresponded to either a free spin {V < Vc) or a Kondo 
temperature which is smaller than the smallest value we 
can reach with our numerical code. In either case, we set 
q — 0, effectively decoupling the free moment from the 
rest of the lattice. Yet, we were still able to span several 
decades of energy scales. 

In the perturbative treatment, the solution of each im- 
purity problem is found by solving a set of two non-linear 
equations for n and /i, which is defined by Eqs. l(3(Hl and 
H35II . As in the slave boson treatment we used the Powell 
hybrid method. The calculation of the second order cor- 
rection for the self-energy involves Fourier transforms as, 
according to Eq. l(3T]l , it has a simpler form in imaginary 
time rather than in frequency space. For this, we used 
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the Fast Fourier Transform algorithmi^i 
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